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Abstract 

Following Ben-Artzi and LeFloch, we consider nonlinear hyperbolic 
conservation laws posed on a Riemannian manifold, and we establish an 
I/^-error estimate for a class of finite volume schemes allowing for the 



0^ 
(N 

(-H I approximation of entropy solutions to the initial value problem. The error 

-i—i . in the L^ norm is of order ^^' * at most, where h represents the maximal 

"3 ' diameter of elements in the family of geodesic triangulations. The proof 

relies on a suitable generalization of Cockburn, Coquel, and LeFloch's 
theory which was originally developed in the Euclidian setting. We extent 
the arguments to curved manifolds, by taking into account the effects to 
s^ , the geometry and overcoming several new technical difficulties. 

o 

v^ ! 1 Introduction and background 

^. . 

f«^ , 1.1 Purpose of this paper 

^rs^ ' The matheniatical theory of hyperbolic conservation laws posed on curved mani- 

(^ [ folds M was initiated by Ben-Artzi and LeFloch [3] , and developed together with 

collaborators (TJ [2l |4l [181 [191 [20j . For these equations, a suitable generalization 
of Kruzkov's theory has now been established and provides the existence and 
S^ ' uniqueness of an entropy solution to the initial and boundary value problem 

H I for a large class of hyperbolic conservation laws and manifolds. The conver- 

gence of the finite volume schemes with monotone flux was also established for 
conservation posed on manifolds. 

The purpose of the present paper is to show that the error estimate for finite 
volume methods, due to Cockburn, Coquel, and LeFloch 7\ in the Euclidian 
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setting carries over to curved manifolds. To this end, we will need to revisit 
Kuznetzov's approximation theory [TBI HI] and adapt the technique developed 
in [7] . One technical difficulty addressed here is the adaption of the standard 
"doubling of variables" technique to curved manifolds. We recover that the rate 
of error in the L^ norm is of order /i^'^, where h is the maximal diameter of an 
element of the triangulation of the manifold, as first discovered in [7j . 

Recall that the well-posedness theory for hyperbolic conservation laws posed 
on a compact manifold was established in [3j , while the convergence of monotone 
finite volume schemes was proved in [T]. In both papers, DiPerna's measure- 
valued solutions [TT] were used and can be viewed as a generalization of Kruzkov's 
theory [TJ]. In contrast, in the present paper we rely on Kuznetsov's theory, 
which allows us to bypass DiPerna's notion of measure- valued solutions. Indeed, 
our main result in this paper provides both an error estimate in the L^ norm 
and, as a corollary, the actual convergence of the scheme to the entropy solution; 
this result can be used to establish the existence of this entropy solution. 

For another approach to conservation laws on manifolds we refer to Panov 
[23] and for high-order numerical methods to Rossmanith, Bale, and LeVeque 
[24] and the references therein. Concerning the Euclidian case M = M" we want 
to mention that the work by Cockburn, Coquel, and LeFloch [5] [7] (submitted 
and distributed in 1990 and 1991, respectively) was followed by important de- 
velopments and applications by Kroner [T3] and Eymard, Gallouet, and Herbin 
[13] to various hyperbolic problems including also elliptic equations. In [B], the 
technique of convergence using measure-valued solutions goes back to pioneer- 
ing works by Szepessy [25l |26] and Coquel and LeFloch [8] [9l [10] . Concerning 
the error estimates we also refer to Lucier [21] [22] , as well as to Bouchut and 
Perthame [5j where the Kuznetsov theory is revisited. 

An outline of this paper follows. In the rest of the present section we present 
some background on conservation laws on manifolds and briefly recall the cor- 
responding well-posedness theory. Then in Section 2 we present the class of 
schemes under consideration together with the error estimate. Sections 3 and 
4 contain estimates for various terms arising in the decomposition of the L^ 
distance between the exact and the approximate solutions. The proof of the 
main theorem is given at the beginning of Section 4. 

1.2 Conservation laws on a manifold 

Let (M, g) be a connected, compact, n-dimensional, smooth manifold endowed 
with a smooth metric g, that is, a smooth and non-degenerate 2-covariant tensor 
field: for each x e Af , gx is a scalar product on the tangent space TxM at x. 
For any tangent vectors X,Y € TxM, we use the notation gx{X, Y) = {X, Y)g 

1/2 

and \X\g := {X,X)g . We denote by dg the associated distance function and 
by dvg = dvM the volume measure determined by the metric. Moreover, we 
denote by Vg the Levi-Civita connection associated with g. The divergence 
operator divg of a vector field is defined intrinsically as the trace of its covariant 
derivative. It follows from the Gauss-Green formula that for every smooth vector 



field and any smooth open subset S C M 



diYg f dvM = / (/, n)g dvQs, 
s Jas 

where dS is the boundary of S, n is the outward unit normal along dS, and 
dvds is the induced measure on dS. 

Consider local coordinates (x*) together with the associated basis of tangent 
vectors {ei} = {di} and covectors {e*}. The differential of a function u : M ^ M. 
is the differential form du = [dujie^ — ^e% where the summation convention 
over repeated indices is used. The vector field V^u associated with du is given 
by VgU — i'^guY Oi = 5*-' {du)jei, where {g^^) is the inverse of the matrix 
idij) — {i^i^^j)g)- The covariant derivative of a vector field X is a (1, l)-tensor 
field whose coordinates are denoted by (VgX)'[. The following formula for the 
divergence of a smooth vector field will be useful: 

divg (/(m, x)) = du{duf{u, x)) + ( divg /) {u, x) 

ax ^\g\ 

We will use the following standard notation for function spaces defined on M. 
For p g [1,00] the usual norm of a function h in the Lebesgue space LP{M;g) 
is denoted by ||/i||Lp(M;g) and, when p = 00, we also write ||/i||oo- For any 
/ G Ll^^{M;g) and any open subset N C M we use the notation 

f{y) dvg{y) := \N\-' [ f{y) dvg{y), \N\g := / dvg. 



g- 

N " JN JN 

1.3 Well-posedness theory 

We are interested in the following initial-value problem posed on the manifold 

ut + divg (/(w, •)) = on R+ x M, (1.1) 

u(0, •) = uo on M, (1.2) 

where u : R+ x A/ — > R is the unknown and the flux / = fx{u) — fiu, x) is a 
smooth vector field which is defined for all a; € M and also depends smoothly 
upon the real parameter u. The initial data in p.2p is assumed to be measurable 
and bounded, i.e. uq £ L°°{M). Moreover, / satisfies the following growth 
condition 

max |( diVg/) (u, x)| <C + C" |u|, ugR (1.3) 

for some constants C, C > 0. 



Definition 1.1. A pair {U,F is called an entropy pair if U : M. ^ M. is a 
Lipschitz continuous function and F = F{u, x) is a vector field such that, for 
almost all u &M. and all x G M , 

duF{u, x) = duU{u)duf{u, x). 

If U is also convex, then {U,F) is called a convex entropy pair. 

The most important example of convex entropy pairs is the family of Kruzkov's 
entropies, defined for u, c G R by 

{U{u,c),F^{u,c)) := (\u-c\, sgn(u - c)(/j,(m) - fx{c))) 

} ' . (1-4) 

= \(u\J c-uKc), fx{u V c) - fx{u A c)j , 

where uM c — maxjit, c}, m A c = niinju, c}. 

Definition 1.2. A function u G L°°(R+ x M) is called an entropy solution 
to the initial value problem (jl.ip - (|1.2p . if for every entropy pair {U,F) and all 
smooth functions (j) = (j}{t,x) > compactly supported in [0,oo) x M, 

(U{u)(f>t + {Fx{u),\/g(f>)g) dVgdt 

^*' (1.5) 

+ // Gx{u)(j){t,x)dvgdt+ U{uo{x))(t){0) dvg>0, 

where Gx{u) := (divg Fx){u) - 3„f7(u)(divg fx){u). 

For instance, with Kruzkov's entropies the above definition becomes (for all 

cG M) 

// (u{u,c)<pt + {Fx{u,c),Vg(j))g)dvgdt 

J JM+xM ^ ' 

sgn(w — c) (divg /)(c, a;) (f>dvgdt + \uo — c\ (f>{0) dvg > 0. 



The well-posed theory for the initial value problem (|l.ip - (|1.2p was established 
in Ben-Artzi and LeFloch [3]. 

In the present paper, we are interested in the discretization of the prob- 
lem (|l.ip - (|1.2p in the case that the initial data is bounded and has finite total 
variation 

uoeL°°{M)nBV{M;g). (1.6) 

In particular, it is established in 3 that in the case of bounded initial data, 
the following variant of the maximum principle is established: 

||w(t)||L-(A./) < Co(r,g) + C^(r,g) |JM(s)|Uoo(M), 0<s<t<T, 

where the constants Cq, Cq > depend on T and the metric g. 



Recall the definition of the total variation of a function w : M 



TVg(w) :— sup / w divgcpdvg, 

\\d>\\r^.<lJM 



-_,<1JM 

where (f> describes all C^ vector fields with compact support. We denote by 

BViM-g) = {ue L\M;g) /TY giu) < ex,}, 

the space of all functions with finite total variation on M. It is well-known 
that (provided g is sufficiently smooth) the imbedding BV{M;g) C L^{M;g) is 
compact. 

In fact, an important property of entropy solutions to (|l.ip - (|1.2p is the fol- 
lowing one: u has finite total variation for all times i > if (|1.6p holds and, 
moreover, 

TV<,(u(t)) < C,iT,g) + C[{T,g) TVg(u(s)), < s < i < T, 

where the constants Ci,C[ > depend on T and g; see |3] for details. Of course, 
this implies a control of the flux of the equation 

sup / divg {f{u{t, •), •)) dvg < C TVg{uo). 

However, as noted in [l], this inequality can be derived more directly from the 
conservation laws and one checks that the constant C is independent of both T 
and g and only depend on the largest wave speed arising in the problem. 

2 Statement of the main result 

2.1 Family of geodesic triangulations 

For r > 0, we consider the uniform mesh t„ :~ nr {n = 0,1,2,...) on the 
half- line IR+ . For h > we denote by T'' a triangulation of the given manifold 
M which is made of non-overlapping and non-empty curved polyhedra K C M, 
whose vertices in dK are joined by geodesic faces. We assume that, if two 
distinct elements Ki , K2 E T'* have a non-empty intersection, say /, then either 
/ is a geodesic face of both Ki, K2 or J{"^^(/) — 0, where !K"~^ denotes the 
(n — l)-dimensional Hausdorff measure. 

The boundary dK of K consists of the set of all faces e of K. We denote by 
Ke the unique element distinct from K sharing the face e with K. The outward 
unit normal to an element K at some point a; g e is denoted by iie^Kix) G TxM. 
Finally, \K\ and |e| represent the n- and (n — l)-dimensional Hausdorff measures 
of K and e, respectively. We set 

PK := Yl 1^1 

eedK 



and for each K dT^ the diameter Hk of K is 

hK := sup dg{x,y). 

x,yGK 

We set 

h := sup{/iK : K e T'*}, 

which is assumed to tend to zero along a sequence of geodesic triangulations. 
We also assume that there exist constants 71 , 72 > such that 

^i^h < T < 71/1 (2.1) 

and 

72'\K\<hKPK<l2\K\ (2.2) 

for all K ^ 7^. This condition implies that (as h ^ 0) 

r-^0, /iV-i->0. 

Finally, we set T = t tit for every integer tit- 

2.2 Numerical flux- functions 

As in the Euclidean case, the finite volume method can be introduced by for- 
mally averaging the conservation law ()1.1|) over an element K G T'', applying 
the Gauss-Green formula, and finally discretizing the time derivative with a two- 
point scheme. First, we define a right-continuous, piecewise constant function: 
for n = 0, 1, . . ., 



u'^it, x) = < (i, x) e [i„, i„+i) X M, (2.3) 



where 



and 



•'K — /I uitn,x)dVg{x), 

Ik 



A:=/ uo{x)dvg{x). (2.4) 



K 



Then, in view of (|l.ip we write 



= — / u{t,x) dvg{x) + I divgf{u{t,x),x)dvg{x) 



u"+i 



- + TW\ E {f{u{t,y)),n,,K{y)),dTg{y). 



eedK' 

We introduce flux-functions fe,K : K x M ^ R and write 



The discrete flux are assumed to satisfy the following properties: 



• Consistency property : for w G R, 

fe,K{u, u) = I{f{u, y),ne,K{y))g dTgiy). (2.5) 



• Conservation property : for u, w G K, 

fe,K{u,v)+fe,KA^,u)^0. (2.6) 



• 



Monotonicity property: 



IUk > 0, lu, < 0. (2.7) 



Then, we formulate the finite volume approximation as follows: 



u 



K ■— ^K 



\K\ 



Y, \e\feM^K,<J (n = 0,l,...). (2.8) 



eedK 



For the sake of stability of the numerical method, we impose a CFL stability 
condition: 

T sup ^ Lip(/) < 1, (2.9) 

where Lip(/) is the Lipschitz constant of /. 

2.3 Main theorem 

The main result of the present paper is as follows. 

Theorem 2.1 (Error estimate for the finite volume scheme on manifolds). 
Let u : R+ x Af ^ R he the entropy solution associated with the initial value 
problem ([LT|) -(fL2 l) for an initial data uq £ L°°{M) Ci BV{M;g). Let u'' be the 
approximate solution defined by p.3|) and ()2.8p . Then, for each T > there 
exist constants 

Co = Co(T,.9,||7.o||l-), Ci=Ci(T,5,TVg(uo)), 

C2 = G2{T,g, ||-Uo||L2(Af;g)) 

such that for all t G [0, T] 

\\u''\t)-u{t)\\Ll(^M;g) 

< (Co |Af|, + Ci) /i+(Co |A/|y2 + (Co Ci)'/') |M|y2/,i/2 
+ ((Co C2f'\M\l/' + (Ci 02^') \M\l/^ h'/\ 



The rest of the present paper will be devoted to the proof of this theorem, 
which will follow from a suitable generalization of the arguments introduced 
earlier in Cockburn, Coquel, and LeFloch [7]. 

Remark 2.2. 1. It is sufficient to establish Theorem 12.11 for smoother initial 
data. When uq is measurable and bounded on AI one can then show the exis- 
tence of weak solutions to the initial value problem (|l.ip - (|1.2p as follows. Let 
u[f G L°°(Af) n BV{M) be such that 

lim Mq = uq in L {AI;g). 

Solving the corresponding problem (|l.ip - (|1.2p for the regularized initial data, 
we deduce from Theorem 12.11 that the approximation solutions {u''}/i>o form a 
Cauchy sequence in L^. Moreover, in view of the (discrete) maximum principle 
established later in this paper and for every T > 0, these solutions are uniformly 
bounded in iy°°((0,T) x M). Consequently, there is a function u E L°°, such 
that 

lim u ^ u in the L norm. 

Finally, we note that Definition 11.21 is stable in the L^ norm. 

2. An immediate consequence of the L^-contraction property is an estimate 
of the modulus of continuity in time, that is 

\\u''{t)-u^'{s)\\mKUg)<Ch + C'\t~s\, 

where the constants C > 0, C" > 1 may depend on T as well as the metric g. 

2.4 Discrete entropy inequalities 

We will rely on a discrete version of the entropy inequality formulated by ex- 
pressing u^^ as a convex combination of essentially one-dimensional schemes. 
For each K E 7'^, e € OK, we define 

^V;' := < - T^(/e,x(UK,<J - /e.K«,0) (2.10) 



■K ~ 


\K\^ 


Je,KyUj^, 


UrJ 


'- J 








rpK 

\K\ 


■Ik, 


Ik 


1 
PK 


eedK 


,k{'u 


"k^y 



and set 
where 



(2.11) 



Therefore, in agreement with (|2.8p we find 

K^'-^H 1^1 "'^e'- (2.12) 



''K 



Remark 2.3. The error estimate will be derived from the family of Kruzkov's 
entropies. Recall that any smooth entropy ri(u) can be recovered by the family 
of Kruzkov's entropies, that is 

The result follows for any entropy by a standard regularization argument. More- 
over, if {rj, q) is any convex entropy pair, then 

q.{u)^\ f dlviOF.Au,OdC 

^ ./TCP 



Define the numerical family of Kruzkov's entropy-flux as 

Fe,K{u,V,c) := fe,K{uy C,V\J C) ~ feM{u/\C,vAc). (2-13) 

Given any convex entropy pair {U, F), define the numerical entropy- flux Fe_K{u,v) 
associated to F by 

F.^Kiu, v):=l f dlU[0 Fe.Kiu, V, d^ 

Hence, from the condition of the discrete flux we see that Fe^x(u,w) satisfies 

• For u e M, 

Fe,K{u, u) = r(Fy{u), ne.K{y))g dTg{y). 

• For u, w e M, 

Fe,K{u,v) -)-Fe,/i:^(w,M) = 0. 

These properties are inherited from the corresponding properties for the numer- 
ical family of Kruzkov's entropy-flux. 

We are in a position to derive the discrete entropy inequalities. For each 
K e T'', e e dK and u,v £R, we define 

where 

rpK 

|A| 
Hence from (|2.10p . iJej^ (u^,u^_^) — v!]^l and by definition of Fi^j^, we have 

duiie.K{u,v)>^, dvHe.K{u,v) >0. 

The last inequality is an immediate consequence of the monotonicity oi fe^xiu, v). 
The former follows from this property and the CFL condition. Moreover, we 



observe that 

He,KiuV X,vW X) - HeJi{uA\,V AX) 

= {u\/X-ZUK {fe,K{u V A, W V A) - /e,/f (w V A, u V A))) 

- [uhX-WK{fe,K{uAX,V AX) - fe^K{uAX,uAX))j 

= (uVA-uAA)- zuk {fe,K{uy A, u V A) - Je,K{u A X,v A A)) 
+ WK{fc,Kiu\/ X,uy X) - fe,K{uAX,uAX)) 

= U{u,X) -WK {Fe,K{u,V,X) - Fe^Kiu,U,X)), 

where we have used (j2.13p . Now, smce iJe,i<'(u,w) is an increasing function in 
both variables, we have 

He,K{u,v)VHe^K{X,X)<He,K{uV X,VV X), 
He,K{u,v) AHe^K{X,X)> He,KiuAX,V AX), 

hence, 

H^^KiuV X,v\/ X) - H^,k{uAX,v AX) 

> {H,^K{u,v)y H,^k{X,X)) - {HeM-^.v) AH,^k{X,X)) (2.15) 

= U{H,,k{u,v),X). 
Consequently, from (|2.14p . (|2.15p taking u — u^ and v = u^ , we obtain 



U{nZ') - f^«) + ^(^e,K«,<J - Fe,;,«,0 



<0. 



Therefore, we have proved: 

Lemma 2.4 (Entropy inequalities for the finite volume scheme). Let {U,F) 
be a convex entropy pair. Then, there exists a family of Lipschitz functions 
Fe^K '■ ^^ ^r ^, called numerical entropy-flux associated to F , satisfying the 
following conditions: 

• Consistency property : for u G M, 

F,^K{u,u)^j{Fy{u),n,^K{v))gdTg{y). (2.16) 

• Conservation property : for u, v G R, 

F,^K{u,v)+F,^KAy.u)=Q. (2.17) 



10 



• Discrete entropy inequality: 

U{uZ') - U{ul) + I|£(f,,^(u^,u'^J - Fe,K«,0) < 0. (2.18) 

From (|2.1ip and (|2.18p , we can write the discrete entropy inequality in terms 
of u^^ and u^, that is 

U{uV:e) - U{ul) + I|^(i^,,^«,<J - F,,K{u'k,ul)) < Dl+l (2.19) 

where i5^+i==f/(«^,i)-t/(«^,i). 

To end this section we also recall the discrete maximum principle established 
in Amorim, Ben-Artzi and LeFloch T: for n = 0, 1, . . . , tit, 

max Kl < ((7o(T) + max K|) C'^{T), 
for some constants Co (T) , Cq (T) > 0. 

3 Derivation of the error estimate 

3.1 Fundamental inequality 

From now on it will be convenient to use the notation Qt '■— [0, T] x M . In this 
section we derive a basic approximation inequality on the manifold Af , that is, 
we derive a generalization of the Kuznetsov's approximation inequality for the 
L^ distance 

\\u\T)-u[T)\\l^(M;s)- 

Before proceeding we need introduce some special test-functions and make some 
preliminary observations. 

Let (^ : E ^ M be any C°° function such that supp(p C [0, 1], f > 0, and 
/ (^ = 1. For each t' e R, a:' € M be fixed, each (5, e > and alH £ M, a; S M, 
we define 

Ps{t; t') -.^ - ^(^^) , ^e(x; x') := - ^{ ^ '' ;, ) , (3.1) 

where we use the Riemannian distance. Observe that pa (t;i') = ps{t';t), 'ip^{x;x') 
ipe{x'; x), and 

ps{t;t')dt = 1, / ■ip^(x;x')dvg{x) ^ 1. 

Clearly, ip^ is a Lipschitz function on M with compact support contained 
in the geodesic ball of radius e, hence ^<; G W^'°°{M) and by Rademacher's 



11 



theorem [T2] it is difFerentiable almost everywhere. Moreover, there exists a 
constant C > 0, such that for iL^— a.e. t gM. and 3-C"— a.e. x G M 

\psit;t')\<j, \dtpsit;t')\<^, 

(3.2) 

Let w : M — > R be a locaUy integrable function on M. We may assume only on 
some N C M and to be extended by zero on M\N. Therefore, there exists a 
sequence of smooth functions {u™} defined on M, such that 

lim w™ ^v on L^{M;g). 

rn — *co 

As in Kruzkov ,15j , in the case of manifolds we can establish the following 
approximation result. 

Lemma 3.1. Given a bounded and measurable function v : M —>■ M. and e > 0, 
the function 

V; := / / 'ipe{x;x') \v{x) - v{x')\ dvg{x)dvg{x') 

J M J M 

satisfies 

lim V, = 0. 

Proof. First, consider the case that v is smooth on M . Let x, x' be two points 
on M and, 7 : [0, 1] ^ M be a minimizing geodesic with 7(0) = x and 7(1) = x' . 
Therefore, for some ^ S (0, 1) one can write 

Hx) - v{x')\ = |«(7(1)) ~ v{im\ = l|(" ° 7)(0| 

= I W7(0), ( J)(0)i < \'^Ai{m I J(o| 

< ||Vgw||oorfg(a;,x')- 
Then, applying the inequality above and using (|3.2D . we obtain 

JM JM 

< j ^\\V gv\\^Vg{B^{e)) edvg{x) 

<C\M\\\^gv\\^e, 

where Vg{^Bx{t)) denotes the volume of the geodesic ball of center x and radius 
e, and the constant C > does not depend on e. 

Finally, suppose that v is measurable and bounded on M . Since M is com- 
pact, V is integrable on M and there exists a sequence of smooth functions {v™} 
defined on M converging to v in L^{M;g). We then conclude with a routine 
approximation argument. D 



12 



For each p,p' G Qt, p ■= {t,x), p' :— {t',x'), we consider the special test 
function (f) defined by 

(t){p;p') := ps{t;t')il;^{x;x'). 

Therefore, as S,e ^ 0, the support of cj) is concentrated on the set {p — p'}. For 
convenience, we introduce the following piecewise approximation of the functions 
ps and V'e- For n = 0, 1, . . ., we define p{t;t'), p'{t;t') as 

P5{t;t') -.^ ps{tn+i;t') (t e [t„,t„+i)), 

(3.3) 

and for all K e 7'', we define ijje{x;x'), ip'^{x;x') by the averages of ipe{x;x') 
along each interface, that is 

Mx;x'):^} ^,{y-x')dTg{y) {x e K,x' e M), 

JdK 

(3-4) 
^[{x-x'):=} Mx;y')dTg{y') {xGM,x'eK). 

JdK 

The following estimate will be useful. 

Lemma 3.2. The functions ipe, tpe be defined by p.ip . and (13. 4p . respectively, 
satisfy the estimate 

sup / \'ip^{x]x) — ^pi{x\x)\dvg{x) < C —, 
xeM J M £ 

where C > does not depend on e,h > 0. 

Proof. Given K, let x' , x be two points on Af and K, respectively. Analogously 
to the proof of Lemma |3. 11 we could write 

JdK 

<\Vgi,,(c;x')\h<C-l^, 

where 

\V gtpe{c; x')\ = sup \ygtpe{y;x')\. 

yedK 

Now, we integrate the above inequality on M and obtain 

\M^; x') - M^; ^')\ dvg{x') <-^h Vg{B,{e)) <C-. 

□ 
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Moreover, we define the corresponding approximations (/)'', cf)^ , d^cj), d^ (j), 
Vg(/) and W'i 4> of the exact value, time derivative and covariant derivative of 
the function (f), respectively: 

(j)''{p;p') ■.^ps{t;t')Mx;x'), (j)'"' {p;p') -.^ p's{t:,t') M^l^'), 

d^c^{p;p') := dtps(t;t') M^;x'), dUiP^P') -^ dr ps{t;t') ^^(x;x'), 

and 

V'lcj){p;p') -.^ ps{t;t')V gi;,(x-x'), 

V'^g ^ip;p') := p'sit;t')VgMx;x'). 

Analogously, for convenient we introduce a piecewise constant approximation of 
the exact solution u, that is 

u(t, x) := u{tn, x), [t e [i„, i„+i), X e M). (3.5) 

Remark 3.3. 1. Note that u represents the entropy solution to problem (|l.ip - 
(|1.2p . while m'' denotes the piecewise-constant approximate solution (|2.3p given 
by the schema p.Sp . By definition we have u^ = u^. 

2. The zero-order approximations of the test function 0, that is, (f)^ and (j)^ 
are due to the explicit dependence of the flux function with the spacial variable. 

3. One denotes by dr ^ V' respectively the time derivative and covariant 
derivative with respect to t' and x' variables. 

Next, let us define the approximate entropy dissipation form 
El^iu.u'^) -.^ ff Ql^{u,u'\t',xy,t',x')dvg{x')dt', 



where 

e'l,iu,u\t',x'y,t',x') 

\uit,x)-u''{t',x')\d^<P 



+ ( sgn(u(t, x) - u'^it', x'))(f{u{t, x),x)- f{u'\t\ x'), x), V'^cf^))) dvg{x)dt 



+ 11 sgn{u{t,x)-u''{t',x')) {divgf){u''{t',x'),x) (f)^ dvg{x)dt 

\u{Q,x)-u^{t',x')\(l){Q)dvg{x)+ / \u{T,x)-u^{t',x')\(t){T)dvg{x). 

(3.6) 
Here, the term Qg ^{u,u''';t' ,x') is a measure of the entropy dissipation associ- 
ated with the entropy solution u. Observe that u defined by (|3.5p appears in 
the first term of the right-hand side of p.6p . This is due to the fact that the 
time derivative of u'' needs special treatment, as was observed in |5]. 
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Analogously, reversing the role of u and u^ ^ we define 



ElSu\u):^ // e1^{u\u{t,x)-t,x)dvg{x)dt, 
J Jqt 

where 

ej,(w'',u(i,a;);t,x) 

{\u'^{t',x')-u{t,x)\dU 

+ ( sgniu'^it', x') - u{t, x)){f{u''{t', x'), x') - f{u{t, x), x'), V^'))) dvg{x')dt' 

+ [[ sgn{u''{t',x') ^ u{t,x)) {div'g f){u{t,x),x') (J)'"' dvgix')dt' 

T 

u''{0,x')-u{t,x)\(f){0)dvg{x')+ f \u''{T,x')-u{t,x)\(j){T)dvg{x'). 

M Jm 

Observing that dt'Ps{t;t') ~ —dtps{t;t'), \/'g^p^{x;x') = — VgVc(a;;x') and 
adding the terms Es,g{u, u^) and Es,e{u'^, u), we get the following decomposition: 

Eliu, u^) + El,{u\u) = Rliu, u^) - Sl,{u, u^), (3.7) 

where 



i?^,(u, u^) := / |u'*(T, x') - u{t, x)\ (l){t, x; T, x')dvg{x')dvg{x)dt 

\u^{t', x') - u{T, x)\ (f>{T, x; t' ,x')dvg{x)dvg{x')dt' 



iQt jm 

,h/,l l\ 
I' 
'Qt J M 

\u^{Q,x') - u{t,x)\ 4'{t,x;0,x')dvg{x')dvg{x)dt 

Qt J M 



\u''{r,x') -^(O,^;)! (t){0, x;t, x')dvg{x)dvg{x')dt, 

IM 

(3.8) 
and 

\u{t,x) -~u^{t\x')\ ■0e(x;x') 
- \u{t,x)-u^{t\x')\ ■0^(x;x')) dtps{t;t') dvg{x')dt'dvg{x)dt 
sgn{u{t,x) -u'^{t',x')) 



Qt J JQt 

hi, I l\ 



^{{f{u{t, x), x) - f{u'\t', x'), x)) ps{t- 1') 

- {f{u{t, x), x') - f(u^{t', x'),x')) p'sit; t'), WgM^; ^')), 
+ (div;, /)(u(t, x), x')(p''' - (divg f){u'\t', x'), x)(f>'') dvg{x')dt'dvg{x)dt. 



(3.9) 
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Passing to the limit as S,e ^ 0, we expect that Rg ^{u, u^) converges to 
\u^'{T,x)-u{T,x)\dvg{x)'- I \u^{Q,x)~u{0,x)\dvg{x), 

M J M 

and, if u is replaced by u and the exact differentials in time-space are used, then 
the term Sg^{u,u^) is expected to converge to zero. 

Finally, we obtain the basic approximation inequality which is derived as a 
lower bound for the term Rg ^{u, u^). 

Proposition 3.4 (Basic approximation inequality). The L^ distance between 
the approximate and the exact solution satisfies 



M 



\u''{T,x)~u{T,x)\ dvg{x) 



<C \u''{0,x)-u{0,x)\dvg{x) 
Jm 

+ C{l + TVg{uo)){e + S) 

+ C sup (Sl{u,u^)+El{u,u'^) + El{u\u) 

0<t<T ^ 

where the constant C > may depend on T and also on the metric g, but do 
not depend on h, e, t, and 5. 

Proof. First, we write (|3.8p as 

Rl,{u, U^) = i?i + i?2 + i?3 + i?4, 

with obvious notation. For i?2, we simply observe that i?2 ^ 0. To estimate 
i?i, we consider the following decomposition 

\u^\T,x')-u{t,x)\ 
= \u\T,x')^u{T,x')\ 

- (\u\T, x') - u{T, x')\ ~ \u'\T, x') - u{T, x)\) 

-{\u'^{T,x')-uiT,x)\~\u\T,x')~u{t,x)\) 
> \u''{T,x')-u(T,x')\- \u{T,x)-u(T,x')\ - \u{t,x) - u{T,x)\. 
Then, using this decomposition in the expression i?i, we have 

i?i = / / \u^{T,x')-u{t,x)\(l){t,x;T,x')dvg{x')dvg{x)dt 
Jq Jm 

>i / \u'\T,x')-u{T,x')\dvg{x') 
^ Jm 

- i ((C + C TVgiuo)) e+{C + C TVgK)) S 
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3. Analogously to item 2, we obtain 

R^>-\ I Iw^O,.^') - u{Q,x')\ dvg{x') 

^ J M 

^\C + C TYg{uo)) e+{C + C TVgK)) 5 



and 



Ra>- I ps{t')\u''{t',x')-u{t',x')\dvg{x')dt' 



'1{{C + C TVgiuo)) e+{C + C TV^K)) s) . 
4. Hence adding all these inequalities, we get 

> I \u^{T,x')-u{T,x')\dvg{x')- I \u^{Q,x')-u{id,x')\dvg{x') 



J M Jm 

-2 II ps{t') \u^{t',x') - u{t',x')\ dvg{x')dt' 

- 3((C + C TVgK)) e+{C + C TV^K)) &) . 

5. Finally, by a simple algebraic manipulation we deduce from the above 
inequality that 

\u^{T,x') - u{T,x')\dvg{x') 

M 



<A + 2 psif)\u''{t\x')-u{t',x')\dvg{x')dt\ 

where 

A= I \u''{0,x')-u{0,x')\dvg{x') 

+ 3 ((Ci(r) + C[{T) TVgiuo)) e + {C^{T) + C[{T) TVgK)) 5 
+ 2 (5t («, u'') + E^u, u^) + El, («, u^)) . 
Then, applying the Gronwall's inequality, we get 

\u^'{T,x')~u{T,x')\ dvg{x') <iA. 

M 

u 

3.2 Dealing with the lack of symmetry 

In this subsection we estimate the lack of symmetry in the term Sg ,{u,u'^). 
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Proposition 3.5 (Estimate of Sgjuju'^)). The following inequality holds 
SI{u,u'^)<C{1 + \\uo\\l^)-^\M\ 

where the constant C > may depend on T and the metric g, but do not depend 
on h, €, T , and 5. 



Proof. Step 1. From (13. 9p we write 

with obvious notation. Consider the decomposition Si ^ S'i+ 5": 

J JQt J JQt 

'ip^{x;x') dtps{t;t') dvg{x')dt' dvg(x)dt, 



S'l^ // Ht,x)-u>'{t',x')\dtps{t;t') 

J JQt J JQt 

i ipeix; x') — tp[{x; x') J dvg[x')dt' dvg{x)dt. 
Using (|3.3p . (|3.4p and the definition of u, it follows that 



\S'i\<ll \ii{t,x)-u{t,x)\ I ip,{x;x')dvg{x') I \dtPe{t;t')\dt' dvg{x)dt 



M 



<T— sup ||u(t) - w(t)||Li(M;ff) sup / / tpe{x;x')dvg{x')dT{x) 

" te[0,T] K^T^JdK J M 



<CT^-{C + CTYg{uo)). 
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M J JQ 



On the other hand, to estimate S'l we integrate by parts 

S'i^-\\ [[ sgniuit,x)-u'\t',x'))utit,x)psit;t') 
J Jqt J Jqt 

(ipe{x;x') — ^p'^{x;x')] dvg{x')dt' dvg{x)dt 
\u{T,x)-u'^{t',x')\ps{T-t') 

( V'elx; x) — i^'^{x; x) j dvg{x')dt' dvg{x) 
\u{Q,x)-u^{t',x')\ps{t') 

( V'e(x;x') — '0g(x; x')\ dvg{x')dt' dvg{x) 
— II II \ut{tTx)\ ps{t\t')\ipf{x\x') — il}'^{x]x')\dvg{x')dt' dvg{x)dt 

\u{T, x) - u'^it', x')\ps{T- 1') + |«(0, x) - u'\t\ x')\ psit')) 



M J JQt 



IMJ JQt 

\ipe{x] x') — i^'^{x; x')\ dvg{x')dt' dvg{x) 
Hence, we conclude that 

h 
e 

Here, with some abuse of notation we have written utdvg{x') to denote the 
integration with respect to the measure ut- 



\S'{\<C-^(T{C + CTYg{uo)) + \M\{C + C\\uo\\L^)). 



Step 2. FinaUy, in order to estimate 5*2 we observe that 



^2=// // sgn{u{t,x)-u\t\x')) 

IQt J JQt 



[{{cfx'{fx{u{t,x)) - fx>{u{t,x)),Vgip^)^ + {dW'g f)^,{u{t,x))ip^) PS 

- {{a,,{,Uu^{t',x')) - U{u'^{t',x'),),\7g^^)^ + {dWg f),{u^it',x'j)^p,)ps 

+ {.U^iuit^x)) - U,{u'^it',x')),Vgij,)g{ps - p's) 

- {dW'gf)ju{t,x)){cl>'' - (/)''')) dVg{x')dt'dVg{x)dt, 



where ax'{,fx{u{t,x)) is the paraUel transport of the vector fx{u{t,x)) from the 
point a; to a;'. We use that / is a smooth vector field on M. For x E M fixed, 
u S R fixed and e > sufficiently small, using the Landau notation O(-) we 
write 



{(Jx' {fx(u)) - fx'{u),Vgtp^{x; x'))g + ( divg /) (u, x') tp^{x; x') 
= {{Vgf)Ju)k{x-x'),VgA{x.x'))g + {tvVgf){u,x')^,{x,x') + 0{dl{x,x')) 
^{{Vgf)Ju),\/g{k{x-x')M^;^')))a + 0(e% 
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where k{x; x') is the tangent vector at x of the minimizing geodesic from x to 
x' . Analogously, we have 

{cy^.,UM))-fin\x),^gM^-x'))g + {d^iYgf){u\x)M^-x') 
= {{Vj){u\x'),Vg{k{x-x')A{x-,x')))g + 0{e'). 

Now, denoting the Lipschitz (1, 1)— tensor field I as 

I{t,x-t\x') := {VgF,>){u{t,x),u''{t',x')), 
we have l^sl < IS'al + IS"^'! + 0(e), where 
^2= // // {l{t,x-t\x'),Wg{k{x-x')ij,{x-x'))) Psdvg{x')dt'dvg{x)dt, 



S'i^ // sgn{u-u'^)({!^,{u{t,x))-,U,{u\t\x')),^g^:){~ps^p's) 

J JQt J JQt 

- (dW'g f)^,{Mt,x)){<l>'' - 4'''')) dvg{x')dt'dvg{x)dt. 



The term S2 is estimated as the final step of item 2, we focus on 6*2 term. 
Applying the Gauss-Green formula with respect to the x variable, it follows 
that 

(j(t,a;';i',a;'),Vg(fc(x;x') V'e(a;;x'))) psdvg{x')dt' dvg{x)dt ^ 0. 

Therefore, subtracting the above expressions from ^2, we obtain 

|5^|<Cr(C + CTV<,K))e, 

where, we have used the fact that, due to the compactness of M and the regular- 
ity of the flux function, the function \7 gf{u, x) is uniformly Lipschitz continuous 
for u in a compact set. Combining this result with the estimation of S2, that is 



h 

-I 
e 

we complete the proof of the proposition. D 



\S'^\ < C T (^^{C + C TYgiuo)) + -{C + C\\uo\\L^)\M\gy 



3.3 Entropy production for the exact solution 

We now consider the approximate entropy dissipation associated with the exact 
solution. 

Proposition 3.6 (Estimate of the quantity Eg^{u, u'')). The following inequal- 
ity holds 

El^iu, u'') < C \M\g (1 + hoik-) + C (1 + TYg{uo))) ^, 

where C > may depend on T and the metric g, but do not depend on h, e, t, 
and 5. 
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Proof. 1. For each c G R and in the sense of distributions we have 

Ut{u{t,x),c) + diy gF^{u{t,x),c) +sgn(w(t, x) - c)(divg /)(c, x) < 0. 

Now, we set c — u^{t',x') for each {t',x') e [0, T] x M fixed. Since u is an 
entropy solution to (|l.ip . for n = 0, 1, . . ., we have 

U{uitn+l,x), U^{t', X')) - U(u{tn,x), u'\t' ,x')))ij,{x] x')dVg{x) 

f {F^{u{t,x),u^'{t',x')),\/gM^;x'))gdvg{x)dt 

tn+1 r 

/ sgn{u{t,x)-u'''{t',x')){dWgf){u''{t',x'),x)-)p^{x]x')dvg{x)dt < 0. 

2. Next, we multiply this inequality by ps{t;t') — ps{tn+i',t') and summing 
the first term in time, it follows that 



IT -J- ^ 

y^ / U{u{tn,x),U^)\ps{tn+l]t')- P&{tn,t')ji}e(.X]x')dVg{x) 



n^ — 1 

{F^{u{t, x), u''), Vgi}},{x- x'))ps{t- 1') dvg{x)dt 
IQt 

+ 11 sgn{u{t,x) - u''){dWg f)^{u'') 'tp,{x]x')psit;t') dvg{x)dt 



+ / U{u{T,x),u''-)ps{T;t')ip^{x;x')dvg{x) 

JM 

- f U{{u{0,x),u^)ps{t')Mx;x')dvg{x)<0. 

JM 

By the definition (|3.5p of u we have the identity 

/ U{u{tn, X), u'^) (ps{tn+i;t') - Ps{tn;t')] 1pe{x; x')dVg{x] 
rtn+1 r 

= / / U{u{t,x),u^)dtps{t;t')^pf_{x;x')dvg{x)dt. 

Jt„ JM 

Therefore, Eg^{u,u'^) is bounded above by 

U{u{t,x),u'')dtP5it;t') 

iipe{x; x') — 'ipt{x; x') j dvg{x)dt dvg{x')dt' . 



QtJJQ 



3. By integrating by parts the above equation with respect to t, it follows 
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that 

/:=// // U{u{t,x),u'')dtps{t;t'){iP,-^P,)dvg{x)dtdvg{x')dt' 
J Jqt J JQt 

= -Y. [ II {u{u{t^+i,x),u'')-U{u{tn,x),u'')){'iij,-i;,) 
„^Q JmJJqt ^ ' 

Ps{tn+i;t')dvg{x')dt' dvg{x) 
U{u{T, x), u'')(V'. - A)ps{T; t') dvg{x')dt' dvg{x) 

M J JQt 

U{u{Q, x), u^){i}j, - A)p5{t') dvg{x')dt' dvg{x), 

IM J JQt 

and thus 

il'T J- p P P 

I<^T // \ut\ps{tn+i;t')\il;^{x;x')~'^p^{x;x')\dvg{x')dt' dvg{x) 

„^0 JMJJQt 

+ 1 II {\u{T,x)-u\t',x')\ps{T-t') + \u{Q,x)^u\t\x')\ps{t')) 

J M J JQt ^ ' ' 

|V'e(a;;a;') — il)f\x\x')\ dvg{x')dt' dvg{x) 
h 
e 

D 



<cI^(t{C + C TVg(uo)) + \M\g (C + C \\uo\\l^)). 



4 Entropy production for the approximate solu- 
tions 

It remains to control the approximate entropy dissipation form for the approx- 
imate solution. 

Proposition 4.1 (Estimate of the quantity Eg^{u'^, u)). The following inequal- 
ity holds 

El^{u\u)<C{l + \\uo\\L^){^+T + e)\M\g 

+ C{l + \\uo\\mM-,,))^\M\y^ 

where C > may depend on T and the metric g, hut do not depend on h, t, t , 
and S. 

Once this estimate is established we can complete the proof of the main 
theorem, as follows. 
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Proof of Theorem \2.1\ 1. First, by (|2.ip there exists 71 > 0, such that t < ji h. 
Moreover, without loss of generahty we can take S = e. Therefore, combining 
Propositions I3.4j | 4. II together and denoting 

AoiT) := C (1 + IImoIIl-), ^i(T) := C (l + TVg(uo)) 

A2iT):^C{l+\\uo\\mM;g)), 

we obtain 

where 

X_i := Ao(r) |M|, /i + A,{T) h + A^i^T) |M|i/2 /ji/2, 

Xo := Ao(r) \M\g h + \\u\0) - u(0)|Ui(Af;3), 
Xi:-Ao(r)|A/|g + Ai(T). 
Then, minimizing with respect to e, we obtain 

2. Next, proceeding as in the proof of Lemma [3.11 and by (|2.3p . that is, 
u''(0,a;) = w^f, we have 

ll^?s:-"o||Li(j\/;g) < / / \uo{z) - uo{x)\dvg{z) dvg{x) 

J MJ K 

< l|VgUo||Li(Af:3) h = TVg(Mo) h. 

Consequently, it follows that 

\\U\T) - u{T)\\L^^M;g) 

< (^Ao{T)\M\g + A,{T)') C h 
+ (Ao(r)|A/|i/2 + (Ao(T) Ai(T))'/') C |A/|i/2 h'/^ 

+ ((Ao(r) A2(r))'/'|A/|i/2 + (Ai(r) A2(r))'/') c jMiJ/^ /^V4^ 

which completes the proof of Theorem 12.11 D 

Proof of Proposition\4l\ 1. Fix isT e T'' and e G dK. For {t,x) G [0,T] x M 
fixed, we set c = u(t,x) and u^{t',x') — u^, for {i',x') S [^n'^n+i) ^ -^> 
(n = 0, 1, . . .) and we define 

i^'U^;x'):=JA{x;y')dVg{y'). 

Therefore, by Definition 13.41 and analogously to Lemma [221 it follows that 

mx;x')-^Pi,{x;x')\dvg{x)<C-, 

M £ 
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where the positive constant C does not depend on h,£ > 0. Now, we write the 
local entropy inequality (|2.19p for K and, since it is also valid for K^, we obtain 
respectively 

Pk 






Llfu{ul+^\c)-U{u]„c) 
Pk ^ 

+ t(f,^k{uI, <_„ c) - Fe,i,«, <, c)) < 



We sum the two above inequalities and from (|2.16p and (|2.17p . we obtain 



Pk ^ ' Pk^ ^ ' 

+ T J^^{Fy,{u-K^,c) ~ Fy,iul,c),n,My'))9drgiy') (4.1) 

Pk PKa 

2. We multiply inequality (|4.ip by |e| V'e e ^^"^ s^™ '^^'^'^ ^11 ^ ^ ^-^ ^'^'^ 

e&dK ^^ KG7'' 



kgt'* 



^E / {Fy'iul,c),n,,K{y'))gi'Ux,x')dr,iy') 

r^^^i. J OK 



e£dK ^^ 
KG7'' 



(4.2) 



where we have used 



E ^C/(«t:e,c)^;e= E ^C/(US:|-;,C)<„ 

^^ Pk, ^„ Pk 

E — f^«.,e,c)<e= E — f/«,e,c)^;e, 
ifGT'' KeT'' 
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Since u^ is a convex combination of m^ ^ and the Kruzkov's entropy U is 
convex, we have by Jensen's inequaUty 

E \K\ uiui+') ^: < E 1^1 ^ u(^te) ^^ 

Kerr'' 
Therefore, from (|4.2p we obtain 

E \K\(u{ul^\c)~U{u-^,c))^[ 

■h 

rJ2 f {Fy'iul,c),n,My'))9M^;y')drg{y') 



Ke7'' 



+ ^ E iFy'iul,c),ne,K{y'))g('ipe,eix;x')-'ipe{x;y')] dTgiy') 

eedK-'" 

Kei'^ 
3. Applying Gauss-Green's formula it follows that 

{Fy,{ul, c), n,^K{y'))g A{x; y') dTg{y') 



(4.3) 



dK 



(div'g Fx'){u'}i,c)ilJe{x;x') dvg(x) + / {F^' {u"]^ , c),'V'gtpe{x; x)) g dvg{x) 



a 

K JK 
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Then, from (14.31) we deduce that 



Y, \K\[Uiul+\c)~U{u-^,c))i,i 



{F,,iu''{t',x'),c),V'M^;x'))gdvg{x') 



"^' / {div'gE,,){u\t',x'),c)Mx;x')dvg{x') 

t'^ JM 

< E ^-^ \e\ €,e DZ' + E ^-^\e\U{uZ\c){^:-^:,^) 

eedK P^ eedK ^^ 

Kerr'' Kerr'^ 

By algebraic manipulation, we see that the expression 
E |Al(t/«\c)-[/«,c))v;: 

''"^^ / (^.'K(t',x'),c),\, 

t^ JM 

t 



{Fx'{u^{t' , x'),c), V' gipeix; x))g dvg{x')dt' 



I sgn{u^{t',x') - c) (divg/)(c, x')V'e(2;;a;') dvg{x')dt' 
is bounded above by 

< r^' I sgn{u''{t',x')~c){div'j){u\t',x'),x')Mx;x')dvg{x')dt' 

Jt'^ JM 

Kerr'' Kerr'' 

Now, we multiply this inequality by p'g{t;t') — ps{t;t'^^-^^) and summing with 
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respect to time variable, i.e. i', we obtain that the expression 

riT — l 



J2 I U{u^K.c)(ps{i:t'n+i)~Ps{t-X))^'Ax:x')dvg{x') 
{F,,{u''{t',x'),c),V'Ae{x-x'))g p'sit;t') dvg{x')dt' 



Qt 

sgniu'^it', x') ~ c) ( div; /) (c, x')Mx; x') ~p's{t; t') dvg{x')dt' 



+ I Uiu'^iT, x'),c)ps{t- T) i,[{x- x') dvg{x') 

M 

U{u^{0,x'),c)ps{t) 'ip'^{x;x') dvg{x') 

M 

is bounded above by 

< // sgn {u'\t', x') - c) ( div; /) {u'\t', x'), x')A{x-, x') p's{t; t') dvg{x')dt' 

J JQt 



eedK 



n=0 eEdK 



"T-1 1^ 



E E ^ 1^1 uiuZ'^'^) (^: - <e) mt') 



n=0 eedK ^^ 



/ Ep^ ^/(F,,«,c),ne,A'(y')>i<e(^;^') 



)-Mx;y') p'sit;t')dTg{y')dt'. 






4. Following the lines of proof of Proposition 13.61 we can also derive the 
identity 

Uiul,c)(^ps{t;t'^^^) - ps{t;0^i;'^ix;x')dvgix') 



t' 



^' / U{u''{t',x'),c)dt'Psit]t')^p',{x]x')dvg{x')dt'. 

t' J M 
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Then, the term Eg^{u, u) is bounded above by 

(E E ^-^\e\p'si^:,.Dl+^')dv,{x)dt 



T „=o eedK 



u'\T, x') - u{t, x)\ps{t; T) + \u^{{), x') - u{t, x)\ ps{t) ' 

""^'T ~ , ' (4.4) 

|'0e(a;; x') — ip'^{x; a;') I dvg{x') dvg{x)dt 

+ II I E / l(div;/)«,2;')-(div;/)«,a;)| 

p'g{t;t') ip^{x;x') dvg(x')dt' dvg{x)dt 

+ 11 I E / \Fy'{ul,c)~ay,iFy{u'k,c))\ 

\i^',{x-x')-xP,(x;y')\ p's{t;t') dVg{y')dt'dvg{x)dt. 

5. We write (j4.4p as E1+E2 + E^ + £4^ + E^ with obvious notation. In order 
to estimate Ei we recall that 

and 

<e'-"^e' = -^/ (div;/)«,x')d^^,(a;'). 

Therefore, since U is convex from Aleksandrov's theorem it has second derivative 
a.e. (see [12]), and we can write 

Dte = ^ ( - duU{ul+lc) ^ ( div; /) (u^, x') dvg{x') + 0{t)) . 
Then, we have 
^1^// / E / l(dK/)«'2;')-(div;/)«,x)|p^^^dz;,(x')(ii' 

+ CT\M\g{l + \\uo\\L^)T 
<C'T\M\g{l + \\uo\\L^){r + e), 

where we have used the fact that for every compact K the function Vg ( divg /) 
is uniformly bounded in K x M . Now to estimate £^2, we observe that 

"E E ^|e|C/«\c)(^^-<,)=0. 

KGT'' 
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Moreover, from [T] we recall the uniform bound 

EV^ " I II n+1 n+l|2 ^ II ||2 , r<ll 

n=0 eedK ^^ 

Key'' 

for some the constant C" > 0. 
Hence, we have 

Kerr'' 



riT — l 






n=0 eeOAT 



n=0 eedK ^ 

Kerr'' 
where we have used p.ip . Applying Cauchy-Schwartz's inequality, we obtain 

n=0 eedK ^^ 

Kerr'' 
<CVTi{T \M\,y/^ ^ (ho|li.(M;,) + C[T)) "^ 

<C2(1+ \\U,\\^.^M;,))^\M\\'\ 

The terms i?3 and E/^ are estimated in the same way that we have already done, 
that is 

E:i<C\M\g(l + \\uo\\L^) ^, 

Ei<CT\M\g(l + \\u4L^^ e. 
Finally, we estimate the last term, that is 

^5- // r E l^lS/ \Fy:{ul.c)-Fy{ul,c)\ 

J JQt jo ^gg-h I-"- I JdK 

\'ij)'^{x;x') -■06 (a;; 2/') I Ps{i\t') dTg{y')dt'dvg{x)dt 
<CT\M\g{l + \\uo\\L^)-^, 
where we have used the condition (12.21). D 
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